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Abstract 



The effective potential of the order parameter for confinement is calculated within the Hamiltonian approach by com- 
pactifying one spatial dimension and using a background gauge fixing. Neglecting the ghost and using the perturbative 
"^q'.gluon energy one recovers the Weiss potential. From the full non-perturbative potential calculated within a variational 
i-H approach a critical temperature of the deconfinement phase transition of 269 MeV is found for the gauge group SU(2). 
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Understanding the phase structure of QCD is one of 
the major challenges of particle physics 1,]. Running and 
upcoming high-energy heavy-ion experiments call for a 
deeper understanding of hadronic matter under extreme 
conditions. The central issue is the understanding of 
the deconfinement phase transition from the confining 
hadronic phase with chiral symmetry spontaneously bro- 
ken to the deconfined quark-gluon plasma with chiral sym- 
metry restored. This transition is expected to be driven 
by the gluon dynamics and lattice calculations show that 
confinement is exclusively determined by the strongly in- 
teracting low-energy gluonic modes [2]. Therefore, under- 
standing the deconfinement phase transition requires non- 
perturbative methods. In quenched QCD reliable results 
on the deconfinement phase transition have been obtained 
by means of the lattice approach [3j, which, however, fails 
at large baryon densities due to the notorious fermion sign 
problem. Therefore, alternative non-perturbative meth- 
ods based on the continuum formulation of QCD are re- 
quired. In recent years substantial progress has been 
achieved within continuum approaches to QCD 0, [H, ID, . 
Among these is the variational approach to the Hamilto- 
nian formulation of QCD in Coulomb gauge 0, Q (see 
also Refs. 0, [l(| for related earlier work). In this ap- 
proach the energy density is minimized using Gaussian 
type ansatze for the Yang-Mills vacuum wave functional. 
Within this approach a decent description of the infrared 
sector of Yang-Mills theory was obtained [ll|, [l2 , 13 , 14\. 



Recently, this approach was extended to finite tempera- 
tures by considering the grand canonical ensemble making 
a suitable quasi-particle ansatz for the density operator 
and minimizing the free energy [lil . [lij ]. In this letter 
we present an alternative Hamiltonian approach to finite 
temperature Yang-Mills theory, which does not require an 
ansatz for the density operator. The finite-temperature is 
introduced here by compactifying one spatial dimension. 



2. Order parameter for confinement 

As is well known, Euclidean quantum field theory can 
be extended to finite temperature L _1 by compactifying 
the Euclidean time dimension to an effective length L. At 
temperature L^ 1 the order parameter for confinement is 
the expectation value (P[A_o]) of the Polyakov loop (P 
path-odering) 
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This quantity is related to the free energy of a static (in- 
finitely heavy) quark at spatial position x. In the absence 
of fermions Yang-Mills theory is invariant under gauge 
transformations U(x° 7 x.) £ SU(N) being periodic up to a 
center element zt £ Z(iV) 



U(L,x) = z k U(0,x) . 



(2) 



Under such a gauge transformation the Polyakov loop 
transforms as P [Aq] — ZkP[Ao] and as a consequence 
(P[Aq]) = in the center symmetric confining phase while 
(P[Aq]) ^ in the deconfining phase with center symme- 
try spontaneously broken. In Polyakov gauge, d$A$ = 0, 
and with Aq residing in the Cartan algebra, in the funda- 
mental modular region P[.Ao] is a convex function of Aq 
and by Jensen's inequality (P[Aq}) < P[(Aq)], instead of 
(P[Ao]), one can alternatively use P[(Ao)] or (Aq) as order 
parameter for confinement [3, Note however, that by 
gauge invariance a non-vanishing (Aq) requires the pres- 
ence of an external background field ao- Choosing ao in 
Polyakov gauge and to satisfy (.Ao) = °o the background 
field becomes an order parameter for confinement whose 
value is determined by the minimum of the effective po- 
tential V [{Aq) = ao]. In Polyakov gauge there are still 
residual gauge transformations satisfying Eq. ([2]), which 
transform Ao to Aq = Aq + fik/L, where /!& is a coweight 
satisfying exp(— /i/j) = Zk, and as a consequence of gauge 
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invariance the effective potential of (Aq) 
the periodicity condition 

V[a +ti k /L] = V[a } . 



ao must obey 



(3) 



This potential was first calculated in Ref. [20J in one-loop 
perturbation theory. It was found that T^[ao] is minimal at 
a = 0, so that (P) ~ P [(A ) = 0] = 1 implying that the 
perturbative theory is in the non-confining phase with cen- 
ter symmetry broken. This, of course, is the expected be- 
havior at high temperatures, where perturbation theory is 
reliable. In this letter we calculate non-perturbativley the 
effective potential V[(.Ao)] m the Hamiltonian approach 
and determine from this potential the critical temperature 
of the deconfinement phase transition. 



3. Finite temperature from compactification of a 
spatial dimension 

Clearly the order parameter (P) m P[(A n )] or (Aq) is 
not directly accessible in Weyl gauge A — 0, which is 
assumed in the canonical quantization. However, by 0(4) 
symmetry, all four Euclidean dimensions are equivalent 
and instead of compactifying the time, one can equally 
well introduce the temperature by compactifying one of 
the spatial dimension, say the a^-axis, and consider (A3) 
as order parameter for confinement. This can be seen as 
follows: 

Consider Yang-Mills theory at finite temperature L~ x , 
which is defined by the partition function 



Z{L)=Ttc- lh ^ 



(4) 



Here H(A) is the usual Yang-Mills Hamiltonian defined 
by canonical quantization in Weyl gauge A = 0. The 
partition function ^ can be equivalently represented by 
the Euclidean functional integral, see for example Ref. 21 1 



Z(L)= J Ij2M„(x)e- 



S[A] 



(5) 



where 



x°— pbc 



L/2 



S[A}= J dx° J d 3 xC{A»-x^) 



(0) 



-L/2 



is the Euclidean action and the functional integration is 
performed over temporally periodic fields 



' '4* 



(7) 



which is indicated in Eq. ((SJ by the subscript a; — pbc. 
This boundary condition is absolutely necessary at finite 
L but becomes irrelevant in the zero temperature (L — » 00) 
limit. 



We perform now the following change of variables 

A° -> C 3 



x° -> z 3 



x 1 -> z° A 1 -> C" 
x 2 -> z l A 2 -> C 1 
A 3 -> C 2 . 



x 3 -t z 2 



(8) 



Due to the 0(4) invariance of the Euclidean Lagrangian 
we have 

C(A' l ,a^) = C(C ,l ,z' i ) (9) 



and the partition function ([5]) can be rewritten as 



Z(L) 



■S[C 



(10) 



z 3 —pbc 



where the action is now given by 



L/2 



S[C^}= Jdz°dz 1 dz 2 J dz 3 £(C^,, 



(11) 



-L/2 



and the functional integration runs over fields satisfying 
periodic boundary condition in the z 3 — direction 



C» (2 



\L/2) = C» 



-L/2) 



(12) 



We can now interprete z° as time and z = 2 2 , 2 3 ) 
as space coordinates and perform a usual canonical quan- 
tization in "Weyl gauge" C° = 0, interpreting C = 
(C 1 ,C 2 ,C 3 ) as spatial coordinates of the gauge field, 
which are defined, however, not on 1R 3 but instead on 
R 2 x S 1 . We obtain then the usual Yang-Mills Hamiltonian 
in which, however, the integration over z 3 is restricted to 
the intervall [—■§,■§]• Let us denote this Hamiltonian by 
H(C,L). Obviously H(C,L ->■ 00) = H(C). Reversing 
the steps which lead from Q to © and taking into ac- 
count the irrelevance of the temporal boundary conditions 
in the functional integral for an infinite time-interval we 
obtain from Eq. ([10]) the alternative representation of the 
partition function 

Z(L) = Tre -f^°H(C,L } = Um Trc -TH(C,L) _ (13) 

T— >oo 

Due to the infinite z - (time-) interval T — > 00 only the 
lowest eigenvalue of H(C,L) contributes to the partition 
function Z(L). The calculation of Z(L) is thus reduced to 
solving the Schrodinger equation H(C,L)tp(C) = Eip(C) 
for the vacuum state on the space manifold IR 2 x S' 1 (L), 
where 5' 1 (L) is a circle with circumference L. 

Let us illustrate the equivalence between Eqs. ([5]) and 
(|13p by means of the free scalar field theory in 1 + 1 di- 
mension defined by the (Euclidean) Lagrangian 



C = \ id^f + ^ 



(14) 



Calculating the partition function for this model from the 
functional integral with the temporally periodic bound- 
ary condition <p(L/2) = <f>(—L/2) one finds 

\nZ(L) = -^Tr\n(-d 2 + m 2 ) 

L 

= -\f dx°J dx^Y, J ^(p 2 n+P 2 + m 2 ) , 
o " 

(15) 

where the p n = 2im/L are the usual Matsubara frequen- 
cies. Representing the logarithm by a proper-time inte- 
gral, carrying out the integral over the spatial momentum 
p and using the proper-time representation of the square 
root one obtains 

In Z(L) = - J dx x E Q (L) , (16) 

where 

n 

is identified as the ground state energy (the lowest eigen- 
value of the corresponding Hamiltonian) of the scalar field 
theory (fT4"l) defined, however, on a compact spatial mani- 
fold S l (L). With the substitution @ x 1 -> z° Eq. JTSJ is 
precisely the representation (fT3|) . 

The upshot of the above consideration is that finite tem- 
perature gauge theory can be described in the Hamiltonian 
approach by compactifying a spatial dimension and solv- 
ing the corresponding Schrodinger equation for the vac- 
uum sector. This equivalence holds in fact for any 0(4) 
invariant quantum field theory. 

The above consideration for the partition function can 
be extended to the finite temperature effective potential 
y[(Ao)]. One finds that V[(A ) — a] can be calculated in 
the Hamiltonian approach from V[(A^) — a] with the z 3 - 
axis compactified. Furthermore, as shown in Ref. [22j | . in 
the Hamilton approach the effective potential ^[(^3) = a] 
is given by the energy density in the state minimizing (H) 
for given (A3), 

Below we calculate the effective potential V [(A3) = 03] 
in the Hamiltonian approache exploiting the representa- 
tion (|13p of the partition function. 

4. Hamiltonian approach in background gauge 

In the presence of an external constant background field 
a the Hamiltonian approach turns out to be most conve- 
niently formulated in the background gauge 

d-A = 0, d = d + &, & ab = f acb a c , (18) 

where the hat " " " denotes the adjoint representation. 
This gauge allows for an explicit resolution of Gauss' law, 
which results in the gauge fixed Yang-Mills Hamiltonian 

H = i /"d 3 a;(J A 1 n(x)J A -n(x)+B 2 (x))+i/ c , (19) 



where n^(x) = — i8/5A%.(x) is the "transversal" momen- 
tum operator (d • II = 0) and 

J a = Det f-D • d) , D = d + A (20) 

is the Faddeev- Popov determinant of the gauge (fTS)) . Fur- 
thermore, 

He = y / d 3 xd 3 y JiV a (x)J A F ab (x,y)p b (y) (21) 

arises from the kinetic energy of the "longitudinal" part of 
the momentum operator. Here 

p a = D • n = - (A - a) • n (22) 

is the color charge density of the gluons, which interacts 
through the kernel 

F = (-D • d) 1 (-d • d) (-D • 1 . (23) 

For a vanishing background field a = the gauge (fT8l) 
reduces to the ordinary Coulomb gauge and H (fT9|) be- 
comes the familiar Yang-Mills Hamiltonian in Coulomb 
gauge [H|. 

We are interested here in the energy density in the state 
4> a [A] minimizing (H) a = (ip a \H\ip a ) under the constraint 
(A) a = a. For this purpose we perform a variational cal- 
culation with the trial wave functional 

rl> a [A] = J A 1/2 4>[A - a] , $[A] = Me Hl^A] ^ (24) 

which already fulfills the constraint (A) — a. For a = 
this ansatz reduces to the trial wave functional used in 
Coulomb gauge However, due to the presence of the 
colored background field the variational kernel w(p) is now 
a non-trivial color matrix. Proceeding as in the variational 
approach in Coulomb gauge Q, from (H) a — > min one 
derives a set of coupled equations for the gluon and ghost 
propagators 

V = (AA)o = \w-\ G=-^((D + a)d)~^ . (25) 

Using the same approximation as in Ref. [l6[ in Coulomb 
gauge, i.e. restricting to two loops in the energy, while 
neglecting He (12T1) and also the tadpole arising from the 
non-Abelian part of the magnetic energy, one finds from 
the minimization of (H) a the gap equation 

u 2 = -A-d + x 2 (26) 

wher^ll 

, , 1 / S 2 \n J\A + a}\ 1 , , s 
^ - "2 ( SA(1)SA(2) ) - 2 * WDGToQ)] 

(27) 



1 We use here the compact notation A(l) = A"* (xi). For Lorentz 
scalars like the ghost, the index "1" stands for the color index ai 
and the spatial position xi. Repeated indices are summed/integrated 
over. 



is the ghost loop (referred to as "curvature" ) with Tq and 
r being the bare and full ghost-gluon vertex. The gap 
equation (f2l)f has to be solved together with the Dyson- 
Schwinger equation (DSE) for the ghost propagator 



G 1 = d • d r (l)GT(2)X>(2, 1) . 



(28) 



Due to the presence of the background field these equa- 
tions have a non-trivial color structure. Fortunately, due 
to the choice of the background gauge (fT5)) . the background 
field enters these equations only in form of the covariant 
derivative d = d + a. Choosing the background field in the 
Cartan algebra the above equations can be diagonalized in 
color space. For simplicity, let us consider the gauge group 
SU(2) so that a — aT^ (the extension to SU(7V) is straight- 
forward). The eigenvectors \a = 0, ±1) of T 3 ab = e a3b 
are the spin-1 eigenstates, see e.g. [24j . Since the explicit 
color dependence is only due to the background field a the 
various propagators have to become diagonal in the ba- 
sis which diagonalizes a. Indeed, one can show that the 
above equations (j2"o) and (12"51) can be consistently solved 
for propagators of the form (in momentum space) 

V aT (p) = 8 aT V a (p) , G aT (p) = 6° T G a (p) . (29) 

In addition, one can show that the propagators T> a (p), 
G" T (p) are related to the propagators in Coulomb gauge in 
the absence of the background field, T>(p), G(p), by 



where 



V°{p)=V{p°), G CT (p)=G(p ff ), 



p = p era 



(30) 



(31) 



is the momentum shifted by the background field. The 
relation (f3U)) implies X a (p) = x(p CT ) an d applies also to 
the transversal projector tki(p) = $kl ~PkPi of the Lorentz 
tensors in Coulomb gauge 

T>ki (p) = t M (p) 1 , xu (p) = t M (p)x(p) . (32) 
2w(p) 

By these relations the gap equation (|2l)]) and the ghost 
DSE (|28p reduces to the ones in Coulomb gauge in the 
absence of the background field [l6| 



^ 2 (p)-p 2 + x 2 (p), 



d- 1 (p) = --I d (p), 
9 



(33) 



(34) 



Here d(p) is the ghost form factor, defined by 

d(p) 



G(p) 



3P 



(35) 



and we have replaced the full ghost-gluon vertex T by the 
bare one Tq, which is known to be a good approximation, 
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Figure 1: The gluon energy u)(p) and the curvature x(p) result- 
ing from the full numerical solution of the variational approach in 
Coulomb gauge as described in Ref. I16|l 



see Ref. [25] . Lattice calculation [26| of the gluon propa- 
gator in Coulomb gauge show that the gluon energy can 
be nicely fitted by Gribov's formula [27| 



;(p) = vV + mVp 2 ■ 



(36) 



A full self-consistent solution of the gap equation (1551) and 
the ghost DSE (|34| reveals that w(p) contains in addition 
sub-leading UV-logs, which on the lattice are found to be 
small. Using Gribov's formula (f31)|) for w(p) and solving 
the gap equation (f53"f for x(p) yields 



X (p) = M 2 /\p\ , 



(37) 



which is indeed the correct IR-behavior obtained in a full 
solution [16] of the coupled ghost DSE and gap equation 
show in Fig. Q] but which misses the sub-leading UV-logs. 

5. The effective potential 

As explained in Seed] the constant background field re- 
siding in the Cartan algebra can serve as order parameter 
for confinement when it is directed along a compactified di- 
mension. Choosing a = ae^ and compactifying the 3-axis 
to a circle with circumference L, the shifted momentum 
(j3Tj) becomes 

p a = p ± + {p n - era) e 3 , p n = 2nn/ L , (38) 

where p ± is the projection of p into the 1-2-plane and p n 
is the Matsubara frequency. In the Hamiltonian approach 
the effective potential of the constant background field is 
given by the energy density in the state minimizing (H) a 



under the constraint (A) = a [22j . Using the gap equation 
one finds for the energy density per transversal degree of 
freedom (H) a / (2V) (V is the spatial volume) in the present 
approximation 



s(a, L) 



1 oo 

El £ 

a n — — co * 



(2nf 



(lo(p°)-x(p1) ■ (39) 
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Figure 2: The energy density J39D as a function of x = aL/(2n). 
L _1 was varied from 260 to 280 MeV (from bottom to top). 



By shifting the summation index n one verifies the peri- 
odicity 

e(a + 2n/L,L) = e(a,L), (40) 

which is a necessary property for the effective potential 
of the confinement order parameter by center symmetry, 
cf. Eq. ©. Neglecting x(p) Eq. gives the energy 
of a non-interacting Bose gas with single-particle energy 
w(p). This quasi-particle picture is a consequence of the 
Gaussian ansatz ((24)) for the wave functional. The quasi- 
particle energy w(p) is, however, highly non-perturbative, 
see for example Eq. (f5B|) . The curvature x(p) in Eq. (I3U1) 
arises from the Faddeev-Popov determinant in the kinetic 
part of the Yang-Mills Hamiltonian (fT9l) . 

In certain limiting cases and for < clL/2tt < 1 the 
energy density (j39j) can be calculated analytically. Ne- 
glecting x(p) and assuming the perturbative expression 
for the gluon energy w(p) = |p| one finds from the 
Weiss potential originally obtained in [2(| 



4tt 2 (aL\- (aL 
euv (a,i) = - F (-j (--1 



(41) 



Neglecting x(p) and using the infrared expression for the 
gluon energy w(p) = A/ 2 /|p| (see Eq. one obtains 



ei R ,(a,i) = 2 



M 2 
L 2 



aL 
2^ 



aL 
2^ 



(42) 



This expression drastically differs from the Weiss potential 
(|41l) : While euv(a,L) is minimal for a = 0, the minimum 
of eiR,(a,L) occurs at a = ir/L corresponding to a center 
symmetric ground state. Accordingly e\jy(a,L) yields for 
the Polyakov loop (P) = 1 while eiR.(a, L) yields (P) = 0. 

Obviously, the deconfinement phase transition is related 
to a change of the effective potential from its infrared be- 
havior ei R (a, L) (j42j) to its UV-behavior e uv (a, L) (jITj) . 
To illustrate this let us approximate the gluon energy uj(p) 
(HO) by 



W (p)«|p|+M 2 /|p| 



(43) 




Figure 3: The Polykav loop (P[a]} evaluated at the minimum a = 
a m i n of the full effective potential shown in Fig. [2] as a function of 
T/T c . 



both, the IR and UV but deviates from it in the mid- 
momentum regime, which influences the deconfinement 
phase transition. With u>(p) given by Eq. (|43p and with 
x(p) = the energy density (|39[) becomes 



e(a, L) = e m (a, L) + euv(a, L) 



f(x) = x 2 (x - l) 2 + cx(x - 1), c = 



3M 2 L 2 

27T 2 



For small temperatures L , eiB,(a,L) dominates and the 
system is in the confined phase. As L _1 increases the cen- 
ter symmetric minimum at x = 1/2 eventually turns into 
a maximum and the system undergoes the deconfinement 
phase transition. In the deconfined phase f(x) has two 
degenerate minima and, starting in the deconfined phase, 
the phase transition occurs when the three roots of f'(x) 
degenerate. This occurs for c = 1/2, i.e. for a critical 
temperature 



T, = L~ 



V3M/ir . 



(45) 



This expression agrees with the Gribov formula (f31>| in 



With the lattice result M = 880 MeV this corresponds to 
a critical temperature of T c ~ 485 MeV, which is much too 
high. This, of course, is not surprising given the approxi- 
mation used to arrive at ([4"5"j) . i.e. neglecting the ghost loop 
x(p) and approximating the Gribov formula by Eq. (|43[) . 
Using the correct Gribov formula (|55|) instead of the ap- 
proximation (|43p only slightly reduces the critical temper- 
ature to T c ~ 432 MeV. It is the neglect of the curvature 
x(p) which pushes the deconfinement phase transition to 
higher temperatures as can be easily seen: From the gap 
equation (l3~3~|) follows that in the deep IR w(p) ap- 
proaches x(p) (S3- Therefore neglecting x(p) in Eq. ([3U1) 
increases the contribution of the confining part ein(a) (|42p 
(relative to that of the deconfining part euv (HUO and thus 
pushes the deconfinement phase transition to higher tem- 
peratures as we will also explicitly see below. 



5 



6. Numerical Results 
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We now turn to a full numerical evaluation of the effec- 
tive potential using for w(p) and x(p) the numerical 
solution of the variational approach in Coulomb gauge ob- 
tained in Ref. by solving the gap equation (j33j) and the 
ghost DSE The results for w(p) and x(p) are shown 

in Fig. [TJ With these results one finds from Eq. (|3l?l) the 
effective potential shown in Fig. [5] From this potential 
one extracts a critical temperature for the deconfinement 
phase transition of T c ~ 269 MeV, which is close to the lat- 
tice predictions of T c = 290 MeV. Let us also mention that 
if one uses for w(p) the Gribov formula (|3"6")) and in accord 
with the gap equation (1331) for x(p) its infrared expression 
(137p one finds a critical temperature of T c ~ 267 MeV, 
which is only slightly smaller than the value T c ~ 269 
MeV obtained above with the full numerical solution for 
to(p) and x(p)- This shows that it is indeed the infrared 
part of the curvature x(p) (neglected in Eq. (pPT| . but fully 
included in Eq. (|3"9"]) and Fig. 0), which is crucial for the 
critical temperature. In view of the ghost dominance in 
the IR this is not surprising. Fig. [3] shows the Polyakov 
loop P[a] calculated from the minimum a m ; n of the poten- 
tial (131)1) shown in Fig. [5J At the phase-transition P[a m i n ] 
rapidly changes from P = to P = 1. 

The value T c = 269 MeV obtained above from the full ef- 
fective potential, Fig. [51 is also close to the range of critical 
temperatures T c = 275 . . . 290 MeV obtained in Ref. [l6| 
from the grand canonical ensemble of Yang-Mills theory 
in Coulomb gauge. It is however not surprising that the 
critical temperatures found in Ref. [16| differ somewhat 



from the value obtained in the pres ent paper. The reason 
is that in the approach of Ref. 16] an additional approx- 
imation is made by using a singe particle ansatz for the 
density matrix. Such an approximation is not necessary 
in the present approach. In this respect the present ap- 
proach is superior over the variational treatment of the 
grand canonical ensemble given in Ref. [lij . 

In Ref. [li| the Polyakov loop potential was calculated 
from a functional renormalization group flow equation ap- 
proach using the Landau gauge ghost and gluon propaga- 
tors as input. For the gauge group SU(2) a critical tem- 
perature of 266 MeV was obtained, which compares well 
with our result of 269 MeV. 

In the present approach the deconfinement phase transi- 
tion is entirely determined by the zero-temperature prop- 
agators, which are defined as vacuum expectation values. 
Consequently, the finite-temperature behavior of the the- 
ory and, in particular, the dynamics of the deconfinement 
phase transition must be fully encoded in the vacuum wave 
functional, as should be clear from the considerations of 
Sec. [3] The results obtained above are encouraging for an 
extension of the present approach to full QCD at finite 
temperature and baryon density. 
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